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ABSTRACT 

Unbiased  estimators  of  the  number  of  individuals  in  a  lot 
possessing  various  patterns  of  types  of  defects  are  constructed. 
Explicit  formulas  are  given  for  cases  of  two  and  three  types  of 
defect.  Application  of  the  formulas  requires  knowledge  of  the 
probabilities  of  various  kinds  of  errors  in  the  inspection  process. 
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1.  Introduction 

In  [1]  [2] ,  we  introduced  the  term  defeat  pattern ,  to  reflect  the  incidence 
of  defects  of  different  types  in  the  same  individual.  If  there  are  k  types  of 
defect,  in  all,  under  consideration,  then  an  individual's  defect  pattern  is 
represented  by  a  binary  (0,1)  sequence  of  k  digits,  with  the  j-th  digit  equal 
to  0,  or  1  according  as  the  j-th  type  of  defect  is  not,  or  is,  present  in  the 

V 

individual.  There  aTe  2  possible  different  defect  patterns. 

We  will  denote  the  defect  pattern  (g^^,.. ,,gk)  by  (g),  and  the  number 
of  individuals,  in  a  lot  of  size  N,  possessing  defect  pattern  (g),  by  D  . 
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where  .  =  £  I  •**!  denotes  sunmation  over  all  2  patterns . 

-  gx*0  g2-o  gk«o 

If  a  random  sample  of  size  n  is  chosen  (without  replacement)  from  the  lot, 
and  Zg  of  these  are  found,  on  inspection,  to  have  defect  pattern  (g),  then,  if 
inspection  is  perfect  a  natural,  (and  unbiased)  estimator  of  D  ,  obtained  by 
equating  observed  and  expected  values  of  Z^,  is  NZ^/n.  If  inspection  is 
imperfect,  the  situation  is  not  so  simple,  but  the  same  method  can  be  used, 
provided  the  probabilities  of  correct  and  incorrect  assignments  of  defect 
patterns  are  known.  Since  the  resultant  estimators  -  ft,  say  -  are  obtained 

s 

by  equating  first  moments  of  sanple  and  population  values  of  the  Z  _'s,  we 

£ 

will  call  them  moment  eetimatore . 

m 

We  will  suppose  that 

(i)  if  a  defect  of  type  j  is  present,  the  probability rchkc  jxiwiu.ii  oe  i  • 
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detected  is  pj ; 
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(ii)  if  a  defect  of  type  j  is  not  present,  the  probability  that  it  will  be 
(erroneously)  detected  is  pj; 
and 

(iii)  for  any  defect  pattern,  decisions  in  regard  to  presence  or  absence  of 

each  type  (or  disjoint  sets  of  types)  of  defect  are  mutually  independent. 

This  last  assumption  may  not  always  be  justified,  but  it  serves  as  a  useful 
point  of  reference. 

With  assumptions  (i)  -  (iii) ,  the  probability  that  an  individual  with  defect 
pattern  (f)  is  classified,  as  a  result  of  inspection,  as  having  pattern  (g)  is 


n  W  Jpt  3  3  (l-pj 

j-i  J  J  J 


(l-gjf.  (1-gJd-fJ 
3  JU-pj)  3  3  ) 


(1) 


We  will  further  assume  that  p^  >  pj  for  all  j .  This  does  not  affect  many 
of  our  formulas,  and  it  is  a  very  reasonable  assumption.  It  merely  means  that 
an  individual  is  more  likely  to  be  classified  as  possessing  defect  type  j  if 
indeed  it  does  possess  this  type  of  defect,  than  if  it  does  not. 


2.  Moment  Estimators 

Our  problem  is  to  estimate  the  D  's  from  the  observed  values  of  the  Z  's. 

I  S 

The  expected  value  of  D  is 

EIZgl  •  UN*1  I(f)D|Pg|f  (2) 

and  so  the  monent  estimators,  Bp  satisfy  the  2^  equations 

nr 
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or,  in  matrix  form 


Z  •  nN_1P  fi 


t- 

where  Z  and  D  are  2  - rowed  vectors  with  elements  (Z^>  tad  (p^>  respectively, 


and  P  is  a  2k  x  2k  matrix  with  pg|f  as  the  element  ir  the  g-th  row  and  f-th 


column.  From  (3)', 


D  -  tof  V*Z 


where  P  A  is  the  left  inverse  of  P.  Since 

E[Z]  »  nN'1?  D 

it  follows  that  E[D]  ■  D,  so  D  is  an  unbiased  estimator  of  p  (E[D  1  ■  D  for 

all  C*).)  8  8 

However  it  is  possible  for  some  of  the  D  *s  to  be  negative  (though  not 

-  2 

all  of  than,  because  I(g)®g  *  N).  Similarly  some  (but  not  all)  of  the  D^'s 
might  exceed  N.  This  latter  possibility  occurs  much  less  frequently. 

For  k»l  (a  single  type  of  defect)  we  have,  of  course 


D1  -  N  -  D0,  and  E[Z1]  -  N^ntDjPj+Q^p'} 

-  N^nOfc’  ♦  D(p^-pp } 


whence 


81  *  N(PrP’i )*1(tt'1Z1-pJ) 


Since  Pj  >  p|,  5^  <  0  if  <  npj 

Dj  >  N  if  Zj  >  npj 

In  practice,  one  would  use  a  modified  estimator  such  as 


(6.1) 

(6.2) 


0  if  <  0 


D{  -  ^  if  0  <  Bx  <  N 


N  if  N  <  D, 


This  estimator  is  not,  in  general,  unbiased. 
Table  1  shows  values  of 

(i)  PrfBj  <  0] 

(ii)  Prfy  >  N] 

(ill)  Bias  (Sj*)  -  E[fi1*]  - 
for  a  few  sets  of  values  of  the  peraeters. 


3.  Distribution  Theory 


Conditional  on  the  numbers  Y  «  {Y^}  of  individuals  with  actual  defect 
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patterns  {(f)}  in  the  sample,  the  joint  distribution  of  the  2  numbers 

{Wg|f}  of  the  Y£  individuals  with  pattern  (f)  which  are  assigned  patterns 

{(g)}  as  a  result  of  inspection  is  multinomial  with  parameters  (Yf’*{Pg|f})  • 

different  (f)  and  (f*)  and  any  (g),  (g*)  (whether  or  not  (g)  =  (g*))  Wg|£ 

*0 

and  Wg*  |  £*  are  mutually  independent  binomial  variables  with  parameters 
(Yf,  Pg|f)»  (Yf* >Pg* |  £*)  respectively.  Since 


For 


zg  ■  kr>  "gif  ■ 

the  conditional  joint  distribution  of  Z,  given  Y,  can  be  represented  symboli¬ 
cally 

Z|Y^*  Multinomial  Ckrtp.ixtt  (9) 

(f)  i  S'i 

*0 

where  *  stands  for  "convolution”  and  ”  means  ”is  distributed  as". 

To  obtain  the  unconditional  distribution  of  Z.  this  conditional  distri- 
bution  has  to  be  conpounded  with  respect  to  Y,  which  has  a  nultivariate 
hypergeometric  distribution  with'  parameters  n;  D;  N,  so  that 


(8) 


Pr[Y  -  jr] 


fN)'1 

N 

UJ  “(f) 

mm* 

(0  <  yf  <  Df  ;I(f)yf  ■  n)  (10) 


From  (8)  and  (9),  moments  and  product -moments  of  the  Zg's  can  be 
evaluated  (see  [1],  equations  (17)). 


4.  Variances  of  the  Estimators 

From  (4)  the  variance-covariance  matrix  of  D  is 

Var(C)  -  (Nn'1)V1Var(Z)  (P_1)T  (11) 

(triiere  the  superfix  (T)  denotes  'transpose'). 


*  7?  .>  o- . '  «,*  it"  «_■  n u..  ^ >■  v;  u-  v.  y- 
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From  [1]  (equation  (17)) 


var(Zg)  «  rpg-n{(n-l)p^  +  (N-n)p^KN-T)’1 


.-1 


(12.1) 

(12.2) 


COV(Zg’Zg*)  *  -nUn-lfrgPg**  (N-n)PgPg*}(N-l) 

*"•  pg  m  N‘^cf)Df?if!  ^  ■  N'^(f)Yiif 

*n4  Vs*  ■  • 

5.  Special  Case:  Uniform  Accuracy  of  Inspection 


If  pj  ■  p  and  pj  ■  p’  for  all  j»l,...,k  -  that  is,  the  accuracy  of  inspec¬ 
tion  is  the  sane  for  all  types  of  defect  -  then 

Pg|f  -  pSV10Cl-p)5oiU-p')So° 

where  sa^  -  number  of  types  (j)  of  defect  fbr  which  fj  ■  b,g^  ■  a. 

For  k»2,  we  have 
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(13.2) 
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Hence 

Dm  •  (Nn'1)Cp-p')"2{p2Z00-p(l-p)tZ01*210)*a-p)2Z11) 

»oi  ■  o»>‘1)(p-p,)'2{-pp'z00n>a-p,)z01*a-p)p'z10-(i-p)(i-p')z11} 

810  ■  o*''1)(p-p,)'2(-pp'z00*(i-p)p'z01*p(i-p’)z10-a-p)a-p')z11} 

Du  ‘  (fta‘1)Cp-p,)"2<pl2z00-p'O-p’HZ01*zl0)*Ci-p')2z11) 

For  exanple,  if  p  ■  0.90  and  p'  ■  0.95  then 

500  -  tbTl{1.121Z00  -  0.125(Z01+Z10)  ♦  0.014Zn) 

Dq1  «  Nn"1C-0.062Z00+  1.183Z01  +  0.007Z1()-  0.131ZU) 

510  -  Nn‘1C-0.062Z00+  0.007Z1Q  +  1.183Z1Q  r0.131Zu) 

®11  *  Nn"1(0.003Z00  -  0.066(Z01+Z10)  +  1.249^} 

As  is  to  be  expected  the  major  coefficient  in  the  expression  for  is  that 
of  Z^,  but  the  contributions  of  the  other  terms  is  not  negligible. 


\u 


-pp’(i-p')  -(l-p)p’  -PP'(I-P’)  d-p)p'd*p')  p(i-p’)  d-p)p’d-p’)  -(l-p)Cl-p’) 

-(i-p)p'2  -pp'(i-p')  -pp'(i-p’)  (i-p)p'(i-p')  (l-p)p'a-p')  pd-p')2  -a-pja-p*) 

p'2d-p')  p'2(i-p*)  p'2a-p')  -p'(i-p')2  -p’a-p’)2  -p'd-p’)2  d-p')3 


For  example 


60U  •  0»>-1)(p-p'r3(pp'2z000-pp,(1-p')(20ol*zoio5-(1-',>p'2zioo 
n>a-p,)2z0U*a-p)p’  a-p')(z101*zu0)-a-p)  a-p')2zm) 

In  these  circumstances,  if  also  D£  ■  2~^N  for  all  (f),  then 
•  Pg  “  2'k(p4p')mC2-p-p')k’m 


(16) 


(17.1) 


where  m(=  m(g))  is  the  number  of  l's  in  (g)  ; 

p*  -  2'kCp2-tp'2)®{l-p)2+(l-p,)2}k"m  (17.2) 

X 

and  p^pg*  -  2  k (p2+p ' 2)m^1  {p  (1  -p ' )  ♦  (1  -p)p ’ }I"01+ini0{ (i -p) 2+ (1 -p ' ) 2 j™00  (17.3) 

where 

ma^(=  ma^(g,g*))  is  the  number  of  types  (j)  of  defect  with  gj  ■  a,  gj  «  b. 

(«11  ♦  "10  *  "01  *  "oo  *  k)- 


Ihe  further  specialization  p'  -  1-p,  corresponding  to  a  constant  probabil¬ 
ity  (1-p)  of  error,  leads  to  further  sinplification.  In  particular. 


pg  -  2‘k  i  ^  -  2*k{p2+(l-p)2}k  for  all  (g) 

p^:  -  2'k{P2*a-p)2>^^{*a-p)>"0l',"1# 


as.i) 

(18.2) 


Note  that  (hiqq  ♦  m^)  equals  the  lumber  of  j ’s  for  which  g^  »  gj;  and 
p2  ♦  (1-p)2  *  l-2p(l-p).  We  will  use  the  notation 

*00  ♦  »u  -  e;  p2+(l-p)2  -  A;  2p(l-p)  -  1-A  -  B 


>*k «  N-n  _  ,-k  n^l  *k» 
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From  (12)  and  (18) 

var(Zg)  -  n.2^(l  -  grj  #  2  'FT 

covtf^*)  »  -n  2  k(^j  •  2  k+  A®^  ®) 


(19) 

I 

(20.1) 

» 

(20.2) 


As  an  example  consider 


®011  "  Q*1  1^2p'1^  6^pC1_P)  2  (zooo+z101+z11q)  "P2  (*"P)  (zi 1 1  +znm  +zni  n) 


000  *101  niOJ 


♦  p3z0ii-(i-P)3Zioo} 


“ill  *001  *010J 
(21) 


obtained  from  (16)  by  putting  p'  ■  1-p  . 
We  obtain  from  (20)  and  (21) 


1  N2, 


var(D011>  “n  *{ 


-  N-9n+8 


(2P-D 


W=TJ 


(22) 


(Details  are  given  in  the  Appendix.  Note  that 


A°(2p-1)’6  -  \  {l*(l-2p)‘V.) 


-2,3 


-k 

The  proportion  DQ,,/N  (»2  in  this  case)  is  estimated  unbiasedly  by 


■/011 


Squ/N  and  this  estimator  has  variance 


3JH  [{l*(l-2p)'2}3  -  Ogpi  ]  -  ^Jjj[{l.(l-2p)'2)3-l]  if  n  < N. 
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APPENDIX:  Derivation  of  the  Variance  of  Dq^ 

From  (20)  and  (21) 

(fta-y2C2p-l)6v*r(B0U)  -  ^S2  -  -  j}$§[A3S2  ♦ 

2(p2  (1-p)4  (3AB2)  +p4  (1-p)  2  (3AB2) 

-p3  (1-P)  3  (B2+2AB2+4AB2+2B3)  -p3  (1  -p)  3B3 

♦P4(1-P)2(3AB2)+P2(1-P)4(3AB2) 

"P  (1"P)  5  (3A2B)  -p3  (1-p)  (3A2B) }  ] 

"  ^  ^ 

-  6{p(l-p)S+2p3(l-p)3-*p5Cl-p)}A2B 
♦12{p2  (l-p)4-»p4(l-p)  2)AB2-8p3(l-p)  3B3] 

“^2  -S^S)  S2-5^}(A3S2-3A4B2+3A2B4-B6) 

vrtiere  S2  ■  sum  of  squares  of  coefficients  of  Z's  in  (Nn”1)”1(2p-1)3D011 

-  3p2 (1-p) 4+3p4 (1-p) 2  +  p6+(l-p)6  -  A3 

-1-1  3  a? 

and  Sj  »  sun  of  coefficients  of  Z's  in  (Nn  )  (2p-l)  D011 

-  3p(l-p)2-3p2(l-p)n>3-a-p)3  -  (2p-l)3 

siting  that  since  A+B  >1, 

A3S2  -  3A4B2  +  3A2B4-B6  »  (A2-B2)3  »  (A-B)3  -  (2p-l)6 
we  obtain  '  > 


,”2,3 


wfyll>  “  H  *  TT  "  l^lf}  “  k  *  ST  [<l+(2p-D'‘}‘,-l  ♦ 


8(n-l 


The  proportion  DQ11/N  is  estiamted  unbiasedly  by  Dq^^/N;  this  estimator  has 
variance 

«l»(2p-l)'2)3  -  1*  W1'  • 
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